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^ ■ A radiative transport model is used to study kinetic equihbration during the early stage 

I— i! of a relativistic heavy ion collision. The parton system is found to be able to overcome 

"5 ! expansion and move toward thermalization via parton collisions. Scaling behaviors show 

I up in both the pressure anisotropy and the energy density evolutions. In particular, the 

^ ' pressure anisotropy evolution shows an approximate as scaling when radiative processes 

, are included. It approaches an asymptotic time evolution on a time scale of 1 to 2 fm/c. 

The energy density evolution shows an asymptotic time evolution that decreases slower 

^ \ than the ideal hydro evolution. These observations indicate that partial thermalization 

O ! can be achieved and viscosity is important for the evolution during the early longitudinal 

^ ' expansion stage of a relativistic heavy ion collision. 

o ■ 

q{ '. 1. INTRODUCTION 

o ; 

0^ . Ideal hydrodynamics is very successful in describing experimental data at the Relativis- 
^ ! tic Heavy Ion Collider. Ideal hydrodynamics requires local thermal equilibrium. The ideal 
• ^ I hydrodynamics equations are valid for systems with at least local isotropy. It is very in- 
^ ' teresting to study the equilibration process to get a more detailed understanding of when 
■ and how much thermal equilibration can be achieved in relativistic heavy ion collisions [ 
dl [21 El in El E] • The isotropization process has been studied with the AMPT model 
In particular, the pressure anisotropy, i.e., the longitudinal to transverse pressure ratio, is 
used to measure isotropization in the central cell in central heavy ion collisions. The string 
melting model is seen to give faster isotropization than the default model due to more 
partons in the initial stage. As the parton-parton cross section increases, isotropization 
increases. The pressure anisotropy crosses one at late times when transverse expansion 
sets in. 

The current AMPT has fixed partonic cross sections and only two to two partonic 
processes. In order to go beyond the current AMPT description, we introduce the medium 
dependent cross sections and particle number changing processes. The gg gg {2 2) 
cross section is taken to be the pQCD cross section regulated by a Debye-screening mass. 
It is inversely proportional to the screening mass squared. The screening mass depends 
on the medium locally and is related to the density of particles. As the density increases, 
the screening mass increases, and the cross section decreases. The gg — > ggg (2 3) 
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cross section is taken to be one-half of the 2^2 cross section, consistent with a more 
sophisticated calculation from Xu and Greiner. The 3^2 reaction integral is determined 
from detailed balance to ensure correct chemical equilibration. Isotropic cross sections 
will be used and should be reasonable for the dense central region of a collision. 

With this radiative transport, chemical and kinetic equilibrations can be studied. The 
evolutions of collision rates are shown in Fig. [1] beginning with 2000 particles in a box of 
5x5x5 fm^ with a temperature of 1 GeV and coupling constant = 0.4. The initial 
2^3 rate is large than 3 — 2 rate and there will be particle production. As time goes on, 
these two rates converge toward the rate in chemical equilibrium. It is interesting to see 
whether kinetic equilibrium is maintained throughout the chemical equilibration process. 
Fig. [2] shows that the energy distribution deviates from an exponential distribution at 
early times. In particular, at 0.3 fm/c, low energy part is closer to the final equilibrium 
value. The system relaxes toward the final thermodynamics equilibrium and at time 3 
fm/c, it is very close to the expected distribution with chemical and kinetic equilibrium. 
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Figure 1. Evolutions of collision rates per 
unit volume. Arrows indicate initial and 
chemical equilibrium rates. 



Figure 2. Energy distributions at different 
times. Dashed lines are initial and chemical 
equilibrium distributions. 



2. RESULTS 



Kinetic equilibration during the early stage of a heavy ion collision can be studied 
with the radiative transport model [ [5]. We will start with either initial isotropic or 
initial transverse distribution. The initial particle distributions are all exponential. Two 
initial energy densities eo = 38, 77 GeV/fm^ will be used. The former corresponds to an 
initial temperature of Tq = 0.5 GeV for the isotropic initial conditions, while the latter 
corresponds to 1 GeV. Two values are used for the strong interaction coupling constant 
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as = 0.3, 0.6. Fig. [3] shows the pressure anisotropy ratio evolutions. With isotropic initial 
conditions, Pl/Pt first decreases due to the infiuence of longitudinal expansion then 
increases as thermalization wins over expansion. In contrast to the case with only 2^2 
collisions, instead of the a^To scaling, there is an approximate scaling. Furthermore, the 
transverse initial condition case has the same asymptotic evolution as the isotropic initial 
condition case. In other words, the memory of the initial pressure anisotropy is lost after 
some relaxation time. As the initial energy density increases or as the coupling constant 
increases, the pressure anisotropy goes closer to one and there is more thermalization 
and more memory loss of initial conditions. The energy density evolutions are shown in 
Fig. m This is for the high initial energy density case, similar results are obtained for the 
low initial energy density case. Fig. H] shows that the energy density evolution at very 
early times is determined by the initial pressure anisotropy and the late time evolution 
is determined by final state interactions. In particular, we observe that the evolutions 
are bounded by 1/r for the free streaming case and when as = 0.6. This falls a 

little short of the ideal hydro behavior and shows that viscosity is important in 

the longitudinal expansion stage. 
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Figure 3. Press anisotropy evolution. Figure 4. Energy density evolution. 
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